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ABSTRACT 

Some of the mass that feeds the growth of a massive black hole (BH) in a galactic center is supplied 
by tidal disruption of stars that approach it on unbound, low angular momentum orbits. For each star 
that is disrupted, others narrowly escape after being subjected to extreme tidal distortion, spin-up, 
mixing and mass-loss, which may affect their evolution and appearance. We show that it is likely that 
a significant fraction of the stars around massive BHs in galactic centers have undergone such extreme 
tidal interactions and survived subsequent total disruption, either by being deflected off their orbit or 
by missing the BH due to its Brownian motion. We discuss possible long-term observable consequences 
of this process, which may be relevant for understanding the nature of stars in galactic centers, and may 
provide a signature of the existence of massive BHs there. 



Subject headings: black hole physics — galaxies: nuclei — Galaxy: center 
dynamics — stars: kinematics — stars: rotation 
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1. INTRODUCTION 

A massive black hole in a galactic nucleus grows by ac- 
creting matter from its surroundings, some of it in the form 
of stars that approach it on low angular momentum orbits 
(loss-cone orbits). When the BH mass to is sufficiently 
small so that its tidal radius r t oc to 1 / 3 is larger than its 
Schwarzschild radius r$ oc to (to < 10 s Mq for the dis- 
ruption of solar type stars), the star is tidally disrupted 
before crossing the event horizon. The accretion of stellar 
debris from such events may give rise to observable "tidal 
flares" (Frank & Rees 1976). Various authors estimated 
the rates, timescales, luminosities and spectra of the flares 
(Ulmer, Paczyhski & Goodman 1998; Ulmer 1999; Syer 
& Ulmer 1999; Magorrian & Tremaine 1999; Ayal, Livio 
& Piran 2000). There is to date only marginal evidence 
for the detection of tidal flares (e.g. Renzini et al. 1995; 
Komossa & Bade 1999; Komossa & Greiner 1999). 

In this letter we propose another tidal process whereby 
the BH may reveal its presence. We consider the fate of 
stars that narrowly escape tidal disruption by the central 
BH after being subjected to extreme tidal distortion, spin- 
up, mixing and mass-loss during periastron passage. The 
total mass in such stars is comparable to that in stars that 
are disrupted, and therefore also to (ft/ fm) m i where ft is 
the mass fraction of the BH that originates in tidally dis- 
rupted stars, and f m is the fraction of a disrupted star's 
mass that is ultimately accreted (f m < 0.5, e.g. Ayal 
et al. 2000). The ratio ft/ fm is significant for a low 
mass BH (to < 10 7 M Q ) that accretes from a low-density 
galactic nuclear core, where mass loss from stellar colli- 
sions is small. Murphy, Colin & Durisen (1991) found 
ft/ fm ~0.15 in Fokker-Planck simulations of the growth 
of a low mass BH in a galactic nucleus (their model 4B), 
assuming polytropic stellar models for estimating the col- 
lisional mass loss. More recently, Freitag & Benz (2001, in 
preparation) found ft/ fm ^0.25 for the same galactic nu- 



cleus model using Monte-Carlo simulations with more real- 
istic stellar models, which are more centrally concentrated 
than polytropes, and therefore lose less mass in collisions. 
They also find ft/ fm as high as ~ 0.65 for low-density 
nuclei where the stellar mass function is weighted toward 
low-mass stars. 

Dynamical analyses of the scattering of stars into loss- 
cone orbits (Lightman & Shapiro 1977; see Magorrian & 
Tremaine 1999 for numeric examples) show that tidally 
disrupted stars in galactic nuclei are typically on slightly 
unbound orbits relative to the BH and that they are pre- 
dominantly scattered into the loss-cone from orbits at the 
radius of influence of the BH, = Gm/ er 2 , where a is 
the velocity dispersion far from the BH. The stellar mass 
that is enclosed within is comparable to m. The scat- 
tering timescale is shorter than the dynamical timescale, 
and so the stars are scattered in and out of the loss-cone 
several times during one orbital period. Following the first 
close passage near the BH the stars are on very eccentric 
orbits with apoastron < 2r/j, and so there is a consider- 
able chance that they will be scattered again out of the 
loss-cone before their next close passage, and thus avoid 
eventual orbital decay and disruption. 

The survival probability is further increased by the 
Brownian motion of the BH relative to the dynamical 
center of the stellar system. In a system composed of 
stars with a typical mass M and radius R and a low mass 
BH, which evolve in an initially constant density core of 
radius r c ~ r^, the amplitude of the Brownian fluctu- 
ations is much larger than the tidal radius, (Ar) /r t ~ 

(r c /R) (M/m) 5/6 > 1 (Bahcall & Wolf 1976). The Brow- 
nian motion proceeds on the dynamical timescale of the 
core, which is comparable to the orbital period of the 
tidally disturbed stars. While these results strictly ap- 
ply only if the BH is embedded in an isothermal system, 
which is not the case in the GC, they are expected to hold 
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generally to within an order of magnitude. The orbits of 
the tidally scattered stars take them outside of r^, where 
they are no longer affected by the relative shift between the 
BH and the stellar mass. Therefore, on re-entry into the 
volume of influence, their orbit will not bring them to the 
same peri-distance from the BH. Both the random motion 
of the BH and the scattering off the loss-cone by two-body 
interactions are expected to increase the survival fraction 
f s to a significantly high value (Fig. 1). More detailed 
calculations, which integrate over the orbital distribution, 
are required to confirm these qualitative arguments. 

The cross section for a close passage with a peri-distance 
< r p scales as r p for stars on nearly parabolic orbits (Hills 
1975; Frank 1978). It then follows that the number of stars 
on nearly loss-cone orbits (r t < r p < 2r t ) is comparable 
to the number of disrupted stars, and so the mass frac- 
tion associated with stars in the volume of influence that 
have undergone extreme tidal disturbance, (ft/fm)fs, is 
significant. 

2. MODEL 

The tidal interaction between the star and the BH is de- 
scribed in the reduced mass system where the star, of mass 
M and radius R, is stationary and the BH approaches it on 
an unbound orbit with a peridistance r p . The tilde symbol 
is used to denote quantities expressed in the units where 
G = M = R = 1. In these units SI = 1 is the centrifugal 
break-up angular frequency and E = 1 is the star's bind- 
ing energy up to a factor of order unity. High rotation is 
the longest lasting dynamical effect of a hyperbolic tidal 
interaction that does not end in capture, and so we will 
use the spin-up of the star by the tidal encounter, Af2, 
as a measure of the long-term effects on the star. A gen- 
eral treatment of spin-up by hyperbolic encounters can be 
found in Alexander & Kumar (2001), where it is applied to 
star-star tidal interactions. Here we apply the formalism 
to the star-BH interaction. 

The tidal disruption radius for slightly hyperbolic orbits 
is rt — fh 1 ^ in the limit fh 3> 1. It is convenient to param- 
eterize the tidal interaction in terms of the the penetration 
factor (3 — r t /r p . The linear tidal coupling coefficients Tj 

are functions of the parameter rj = r^j" 1 / y/l + fh , where 
r\ ~ /?~ 3 / 2 in the limit fh ^> 1, and of the eccentricity 
of the orbit e = 1 + 2r p E /fh, where E a is the orbital 
energy of the 2-body reduced mass system. The spin-up 
due to the transfer of angular momentum from the orbit 
by the tidal interaction is related to the energy transfer by 
AE = ItlpAil, where £l p is the relative angular velocity at 
periastron, / is the star's moment of inertia, and rigid body 
rotation is assumed (Goldreich & Nicholson 1989; Kumar 
& Quataert 1998). On parabolic orbits A£l is independent 
of m to leading order in the multipole expansion, and can 
be expressed as (e.g. Press & Teukolsky 1977) 

Ao^ T2(r i /2 y / 2 . (i) 

y/2I 

The linear analysis of the spin-up breaks down as r p de- 
creases and AS1 approaches 1 . At first Afi increases faster 
than predicted by the linear analysis, but then it peaks at 
the onset of mass-loss as the ejecta carry away the angular 
momentum. In the results presented here we do not cor- 
rect for the non-linear effects, and so formally very high 



values of Att should be understood as indicating any com- 
bination of mass- loss, mixing or spin-up effects. The tidal 
coupling coefficients can be calculated numerically for any 
given stellar model. In our analysis we will use two stellar 
models: a, n — 3/2 ideal gas polytrope, and a solar model 
(Christensen-Dalsgaard et al. 1996). 

We demonstrate the tidal scattering effect by consider- 
ing the fate of a star on a circular orbit at an initial radius 
r% = Th that loses angular momentum, but not energy, in 
a 2-body interaction and is deflected into a near loss-cone 
orbit with peri-distance r t < r p <C r t . We assume for sim- 
plicity that the stellar mass within n equals the BH mass, 
that it is negligible within r p , and that its distribution is 
proportional to r~ a (a ^= 2, 3). It then follows from en- 
ergy conservation on orbits in a spherical mass distribution 
that the orbital energy of the reduced mass system at r p 
(i.e. excluding the potential of the stellar cluster) is 
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and the eccentricity is 

e = 1 + 
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3. RESULTS 

We apply our model to the 3 x 10 6 M BH in the Galac- 
tic Center (GC), which is of particular interest because 
present-day observations can already resolve individual 
stars very close to the BH (Genzel et al. 1997; Ghez et 
al. 1998). The BH in the GC is surrounded by a 2-body 
relaxed cusp with a ~ 1.5 and has a radius of influence 
of r h ~ 3.6 pc (Alexander 1999). For the BH in the GC, 
r t > 10rg for a solar type star, and so general relativistic 
corrections can be safely neglected (Laguna et al. 1993). 
Figure 2 shows the cumulative values of the ratio between 
the rate of near loss-cone encounters (r p > r t ) and the to- 
tal rate of disruptive encounters, the orbital period after 
the first periastron passage, and the apoastron distance af- 
ter the first periastron passage, all as functions of A(l. The 
orbits have E Q = 0.04, but are still effectively parabolic, 
with 1 < e < 1 + 10~ 5 . For the solar stellar model, the 
orbital energy remains positive after the first periastron 
passage for all values of AQ < Aflt, = 0.43 (r p /r t > 1.13), 
which corresponds to unbound and tidally disturbed stars 
(defined here as having 1/2 < AQ/A^li, < 1) of total mass 
0.16 (ft/fm) f s tn. For the polytropic stellar model, which 
is less centrally concentrated and therefore more easily per- 
turbed by the tidal interaction, the orbital energy remains 
positive after the first periastron passage for all values of 
Aft < 0.28 (r p /r t > 1.70), which corresponds to unbound 
tidally disturbed stars of total mass 0.19 {ft/fm) fsfn. The 
plots of the orbital period and apoastron radius confirm 
that the qualitative arguments presented in §1 hold for all 
but the closest periastron passages. 

4. DISCUSSION 

We have argued that a significant fraction of the stars 
around a supermassive black hole in a galactic center expe- 
rience considerable spin-up, due to extreme tidal interac- 
tions with the central BH. Unlike the star-star tidal inter- 
actions considered by Alexander & Kumar (2001), which 
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also lead to high rotation, this mechanism does not re- 
quire an extremely high stellar density. We further note 
that the tidally disturbed stars will remain in the galactic 
center even after tidal scattering ceases in the later phase 
of the BH growth when its mass exceeds ~ 10 8 Mq and 
rs > r t . Assuming f s <~0.5, ft/fm ^0.25 and solar type 
stars, and considering only orbits that are unbound to the 
BH after the first periastron passage, we estimate that the 
Galactic Center contains 10 4-5 stars that have survived 
an extreme tidal interaction with the central BH. 

We now examine briefly potential observational signa- 
tures of such interactions. The effects of rotation on stellar 
evolution have been studied extensively by a number of au- 
thors (see e.g. Maeder & Meynet 2000; Heger & Langer 
2000, and references therein). Here we want to concentrate 
in particular on a few observable effects: (f ) The luminos- 
ity, (2) abundance anomalies, and (3) colors. Rotationally 
induced mixing results in higher bolometric luminosities 
for a given mass, due to the increase in the mean molecu- 
lar weight ii, and in the convective core size. Mixing also 
brings the products of hydrogen burning to the stellar sur- 
face. Thus, one can expect enrichment in 4 He, 14 N, 13 C, 
and 26 Al, and depletion in 12 C, 16 0, 15 N. Combining these 
two effects, one of the consequences of high rotation rates 
is stronger abundance anomalies in more luminous stars. 
Interestingly, these were precisely the effects observed by 
Carr, Scllgren & Balachandran (2000), in their abundance 
measurements of the M supergiant IRS 7 in the Galactic 
center. The low C and O and high N abundances indicate 
that IRS 7 has undergone mixing in excess of the values 
predicted by standard models (or observed in other super- 
giants like aOri). In fact, the observed abundances have 



prompted Carr et al. to conclude that "extra mixing in- 
duced by rapid rotation may indeed be the fundamental 
difference between the evolution of massive stars in the 
Galactic center and those elsewhere in the Galaxy." The 
third effect is relevant particularly to stars that would have 
otherwise (with no tidal spin-up) been red supergiants or 
Asymptotic Giant Branch (AGB) stars. The envelopes of 
such stars (particularly the latter) are relatively loosely 
bound to their cores. Consequently, even modest amounts 
of tidal interaction and spin-up lead to a loss of all or a part 
of the envelope (e.g. Livio 1994). The stripped, more com- 
pact star will thus become blue. In fact, this mechanism 
has been proposed to explain the existence of a strange 
group of very blue stars in the globular cluster MI5 (De 
Marchi & Paresce 1996), and the colors of Sk -69° 202, 
the progenitor of SN1987A (e.g. Woosley 1988). Our cal- 
culations therefore predict that the stellar population in 
galactic centers will exhibit a paucity in red giants, and 
will rather be relatively rich in blue stars. Concentrations 
of blue stars around the central BH are found in both the 
Galactic Center (Genzel et al. 1997) and in M31 (Lauer 
et al. 1998). Many of the blue stars are predicted to be 
on elongated ("radial") orbits. Observations of the Sgr A* 
stellar cluster in the Galactic center appear to be consis- 
tent with this prediction (Genzel et al. 2000), although 
the orbital solutions are still quite uncertain (Ghez et al. 
2000). 

We thank Marc Freitag, Andrea Ghez, Tod Lauer and 
John Magorrian for useful discussions. TA is grateful for 
the hospitality of the Aspen Center for Physics, where this 
work was initiated. 
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Fig. 1. — A schematic representation of the change in the disruption probability of a tidally disturbed star as function of the first peri- 
distance r p i when stochastic processes are included. For an isolated 2 body system (i.e. no other stars present) there exists a critical first 
peri-distance r c \ such that a star passing at r v < r c \ will ultimately be tidally disrupted, and conversely, a star passing at r p > r c i will 
survive subsequent close passages. When processes such as scattering and the BH Brownian motion arc included, some of the stars with 
r p < r c \ will survive (hatched area). 
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Fig. 2. — (a) The ratio between the rate of stars that undergo extreme tidal disturbance and those that are tidally disrupted in a model 
of the GC (Mbh = 3 x 10 6 Mq, n = = 3.6pc, a = 1.5). The arrow marks the range of AS7 where the star remains unbound to the BH 
after the first pcriastron passage, (b) The orbital period after the first periastron passage. The arrow marks the orbital period at the radius 
of influence, (c) The apoastron distance after the first periastron passage. The arrow marks the BH radius of influence. 



